Effective quantum equations for the semiclassical 
description of the Hydrogen atom 



Guillermo Chacon-Acosta 

Departamento de Matematicas Aplicadas y Sistemas, Universidad Autonoma 
Metropolitana-Cuajimalpa, Mexico D. F. 01120, Mexico 

E-mail: gchaconScorreo . cua . uam . mx 
Hector H. Hernandez 

Universidad Autonoma de Chihuahua, Facultad de Ingenieria, Nuevo Campus 
Universitario, Chihuahua 31125, Mexico 

E-mail: hhernandezOuach.mx 
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of quantum mechanics. This treatment allows us to describe the quantum state of 
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configuration variables, their moments and quantum dispersions. It also provides a 
semiclassical description of the orbits and the evolution of observables and spreadings 
and their back-reaction on the evolution. 
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1. Introduction 

It is known from general principles [I] that the evolution of the expectation value of the 
position operator in a general quantum state is given by Ehrenfest's theorem, which is 
the quantum analogue of Newton's second law. It can be written as fn'^^^^ = (^F (x)^. 
We notice immediately that although the left hand side corresponds to the classical 
expression, the right hand side does not, instead we have a non-local force term which, 
in principle, can be interpreted statistically. When we expand this term in Taylor series 
we notice that it contains infinite quantum corrections that are expressed as powers of 
the dispersions: 

(x)) = F ((x)) + ^Ax^F" ((x)) + . . . 

The momentous quantum mechanics is based on the same idea applied to any 
operator [2]. It is similar to the low energy effective action method which has been 
largely and successfully used in quantum field theory. There, the actions of interacting 
theories can be seen as quantum corrections to the classical action |3]. Indeed, in some 
circumstances, the effective equations reproduce the results of the effective action even 
better than the well known WKB approximation jl]. 

The effective equations approach for quantum systems has been developed to 
systematically analyze quantum effects through quantum induced corrections to classical 
equations, leading to observable phenomena as deviations from the semi-classical 
behavior [2]. These quantum corrections come from quantum backreaction effects. It 
has been successfully applied to an isotropic and homogeneous model in loop quantum 
cosmology [5], which has been used to predict the evolution of the universe before 
the Big Bang [B]. It has also been successfully used to analyze cosmological isotropic 
models with matter [7] and with positive cosmological constant [8], and also to study 
the effective constriction equations for loop gravity and relativistic systems [H]. 

Effective equations describe the behavior of the expectation values in a definite 
state, replacing the description in terms of the Schrodinger equation by a system of 
infinite coupled equations for both classical variables and quantum fluctuations. These 
equations are particularly suitable for semi-classical states, which may shed some light 
on effects that could be potentially observable. Moreover, in those states, the system of 
infinite equations can be reduced to a finite one by making some consistent truncations, 
as for example in adiabatic approximation. 

For some systems like the harmonic oscillator |2] it is possible to determine the 
exact properties in some state because the second-order moments form a closed system 
and because they decouple from the classical variables. For more complicated systems 
(e.g. anharmonic systems) the equations for both classical and quantum variables are 
highly coupled, so it will require additional approximations for the equations. 

In this work we use the momentous method to obtain the effective description of 
a standard spinless quantum system with a Coulomb potential, i.e. a hydrogen atom. 
In section [2] we describe this version of quantum mechanics with the introduction of 
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the moments or fluctuations as additional variables that encode the quantum degrees 
of freedom. We describe the construction of the quantum Hamiltonian, first for one 
dimensional systems and also for higher dimensions. We also compute the equations of 
motion for the expectation values of basic operators for which the main ingredient is the 
calculation of the Poisson algebra between quantum variables. In section [3] we compute 
the quantum Hamiltonian for this system that corresponds to the classical Kepler 
problem with correction terms. With this Hamiltonian we write the corresponding 
effective equations up to second order in the moments, finding twelve coupled equations: 
three for the classical variables and nine for all the relevant moments. In section |4], we 
solve numerically the system of effective equations and analyze the quantum corrected 
behavior of the system comparing it qualitatively to the classical case. Finally in section 
[5] we summarize and discuss our results. 



2. Effective dynamics of quantum systems: Momentous quantum mechanics 

Effective equations of quantum systems describe the dynamics of expectation values of 
observables, as well as the evolution of its dispersions. These equations allow us to study 
the quantum evolution by analyzing how quantum effects modify classical dynamics. In 
regimes where the fluctuations and dispersions are small with respect to observables 
one can treat the quantum effects perturbatively thus providing the ideal scenario to 
perform numerical analysis in a semiclassical regime [2] . 

Equations of motion in effective theory are derived from an effective Hamiltonian, 
treating expectation values of observables and their associated momenta as classical 
canonical variables. Quantum fluctuations are introduced as a set of new quantum 
dynamical variables defined as follows: 

G'^'' = ((£-^r(p-p)V,, a + h>2 (1) 

where x = (x) and p = (p), and the subscript indicates that the operators inside the 
brackets are Weyl (or completely symmetrical) ordered. This quantum variables describe 
the spreading of the quantum modified evolution from the classical one. 

The momenta are not arbitrary but subject to generalized uncertainty relations 
such as 

^2.0^-0.2 _ > ^. (2) 

Notice that G"^'^ and G^'^ are the standard dispersions Ax"^ and Ap^, and that ^ 
simplifies to the usual uncertainty principle for pure states [Tj. 

Evolution is obtained by evaluating Poisson brackets of variables and momenta 
with the quantum effective Hamiltonian that is defined as the expectation value of the 
standard Hamiltonian operator 

<^>-^« = ^(-^) + Ei|^G- (3) 

a,b ^ 
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where H{x, p) is the classical Hamiltonian. Equations of motion are obtained in the usual 
Hamiltonian formulation / = {/, ifg}. Quantum variables G"'''' are now dynamical, as 
the classical ones {x,p). For general models one obtains an infinitely coupled system 
of equations with infinitely many variables that, however complicated, provides us a 
full description of the system. It is also important to note that the effective equations 
so obtained are state dependent, since the dynamic quantum fluctuations affect the 
behavior of the expectation values. 

To obtain effective equations of motion one compute the Poisson brackets of 
variables with the quantum-corrected Hamiltonian: the following relations, for one 
degree of freedom, are useful 

{x,p} =1, (4) 

{x,G^''} = 0, (5) 

{p, G'^n = 0. (6) 

The Poisson algebra for the moments was originally obtained in [2] but was recently 
reexamined in [8] for its algebraic and numerical implementation 

fj ^.V.G"^^^-"'''^'^-", (7) 
,1. ' 

where the sum runs over odd numbers from n = 1 . . . N, with 1 < < min[a + c,b + 
d,a + b,c + d\, and the coefficient is 

n 




For the quantum harmonic oscillator it was shown that the ground state energy is 
added to the classical Hamiltonian [2] , [10] . It was also seen that the system is solvable 
since the moments are not coupled with the expectation values of classical variables: 
there is no quantum back-reaction. In the case of anharmonic systems, one can use an 
adiabatic approximation that produces effective forces coming from the coupling terms 
between the expectation variables and the moments. This method has also been applied 
to the isotropic and homogeneous model in loop quantum cosmology [5l|6], and to study 
different cosmological models with matter [7], cosmological constant [8], among others. 

For the case of k pairs of canonical degrees of freedom we have the general definition 
for quantum variables: 

^a„6„...,a„...A ^ _ 3,^)^1 (p^ _ p^f. ...(£,_ 3,,)<^.(p^ _ p,f>^)^^^^ (9) 

where + 6^ > 2 and i = 1, . . . , A;. We have GO,o,...,o,o ^ QO,...,a„...,o ^ q_ ^^^^ 
also note that we will always have an even number of indices and each pair (oj, 6j) 
corresponds to the moments of each canonical pairs {xi,pi). 

The generalized quantum-corrected effective Hamiltonian is as follows 

HQ:=y---y — ^ — — r T- G'''''''-''''^'^ (10) 



ai,bi Ofeb, 



ai!6i! ■ ■ ■ Ofclfofc! dxl^dp\^ . . . dxl^'dpl 



Effective quantum equations for the semiclassical Hydrogen atom 



5 



With this Hamiltonian we will obtain the evolution equations for each degree of freedom. 
For this we need the set of Poisson brackets: those among the expectation values are 
known and those between the expectation values and moments vanish. 

The Poisson algebra among the moments is the generalization of equation ([t]); this 
expression was first found in [2] and then corrected in [S] 

k 



i=l 



^,^,Qai,bv,—;ai,bi-^, — \o,k,bkQCi,di;...;Ci-l,di;...;Ck,dk'^ 

+EE E(-i)'(f)""^E..„. 

n s ei,...,ek 

j^n,s,{e} ^a^+ci-ei,;)i+di-ei;...;afc+Cfc-efc,f)fe+dfc-efc 
^ '^{a},{b},{c},{d} ^ 



n = 1, . . . , N , and 



1, ^^(min[ai,di] +min[6i,Q]) < 1, 

'^■{mm[ai, di] + min[6j, q]) - 1, '^■{mm[ai, di] + min[6i, q]) > 1. 

s = 0, . . . ,n ; < Ci < min[aj, d^, s] + min[6j, q, n — s]. 
The /C coefficient is 



j^nMe} _ V- "E..,.n-s TT X ' 9^ J [ 9^ J [ 9^ J [ e. - 9. ^ 

'^{a},{b},{c},{d} - _ „^| 11 / \ / \ ' y^'^) 



s'An — SI 



n — s \ I s 
9i / V 9i 



where max[ei — s, Cj — a,, Cj — di, 0] < gi < min[6j, q, n — s, e,]. 

Henceforth we will consider the case k = 2, corresponding to a two dimensional 
problem, being the case for the Kepler system. 

3. Effective quantum Hamiltonian for Coulomb potential 

In non-relativistic quantum mechanics the hydrogen atom is studied with the 
Schrodinger equation with a central potential; when one considers stationary solutions 
it is possible to separate the angular and the radial parts and thus obtaining the 
corresponding eigenfunctions. It is well known that for the Coulomb potential V{r) = 
— ^, with k = the radial function is the product of the associated Laguerre 

polynomials by a decaying exponential of r and [Ij. Moreover, the corresponding 
energy levels of the Hydrogen atom are = — ^ ^ 2 in terms of the Bohr radius 

The Schrodinger treatment of the Hydrogen atom has long been known and solved 
[IT] . The state of any quantum mechanical system in this picture is encoded in the 
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wave function ?/'(r, t), observables and expectation values are all obtained from the wave 
function considering it as the probability distribution for the system at hand. 

We study now this system using the effective momentous method of quantum 
mechanics exposed in the previous section. We will analyze how the moments evolve in 
time and how they modify the corresponding classical dynamics. 

The classical Hamiltonian for the Hydrogen atom is that of the Kepler problem 



H{r,p) 



2m 2mr2 



k 



(14) 



which, by virtue of the conservation of angular momentum, is a two dimensional 
system with classical polar configuration variables {r,6) and {pr,Pe) are their canonical 
conjugate momenta. 

Using the definition ^ of the moments for our two pairs of canonical variables we 

have 



G 



a,b,cA 



PeY) 



Wcyl 



(15) 



from ( 10 ) we obtain the corresponding quantum effective Hamiltonian 

^2 



I, /^0,2,0,0 r^0fifl,2 
Ho = ^ + - - + — + + 2^ 

a>2 



n 

2m 



r 



2m 



2mr^ 



2mr 



-k 



a+l 



-G 



a,0,0,0 



b>i 



-V. 



mr 



6+2 



ib + l) 



PeG 



fe,0,0,l 



1 



-G 



b,0,0,2 



(16) 



The first three terms correspond to the classical Hamiltonian (14), the fourth term is 
the radial momentum dispersion, the fifth term is the angular momentum dispersion. 
The sum for a > 2 is associated with the fluctuations in the radial component, while 
the sum for 6 > 1 contains the coupled moments for the fluctuations in r and p. 

As can be seen from (16) the system is highly coupled, in contrast with the case of 
the harmonic oscillator where the moments decouple from the expectation values and 
its contribution to the quantum Hamiltonian is constant, rendering the ground state 
energy ^ [2] . The fact that for the central potential the moments can not be decoupled 
may hinder the calculation, however we will perform consistent truncations at different 
orders in quantum variables. We note that, to second order in the moments, the last 
term of the second sum does not appear, that is, in order to consider all the terms 
of the effective Hamiltonian to this order we need to consider at least the third order 
momenta. 

We obtain the equations of motion for this set of variables by computing their 
Poisson brackets with the Hamiltonian (16). For expectation values 

Pr 

m' 

pI k + 

Pr = iPr,HQ^ = :^,--, + —r + J2( 

a>2 



r 



{r,HQ} = 

■ {Pt^Hq} 

{h+l){h + 2) 



mr 



G 



a,0,0,0 



mr-^ 



(17) 
(a + 2) 



2m 



fe>i 



mr 



6+3 



PeG 



6,0,0,1 



4- 



6,0,0,2 
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a>2 



mr 



{pe,HQ} = 0. 



(20) 



Equation (17) is just the usual definition of the momentum associated with r while 
equation (20) states the classical conservation of angular momentum, i.e. po = I = const. 
One can see the strong quantum back-reaction of this semiclassical approach for classical 
variables in equations (18) and (19). 



The equations of motion for moments at any order follow from (11): 

1 



^01,61,02,62 _ 1^01,61,02,62 Hq^ 



2m 



01,61,02,62 ^0,2,0,0 



o>2 >- 



a + 1) 



2mr 



2mr 

r (^01 ,61 ,02, 62 ^o,o,o,ox 

».o+l L ' J 



I _j_ r^oi,6i,02,62 ^0,0,0,2| 



6>1 



mr 



I (^01 ,61 ,02, 62 ^6,0,0,l| _| 1^01,61,02,62 ^6,0,0,2| 

2 



(21) 



Pr 



m 



Pr 



^ _l_ ^ |_ _;^^2,o,o,o 



e = — 

mr 
^1,1,0,0 

^2,0,0,0 
^0,2,0,0 

^0,0,1,1 

(^0,0,2,0 
^0,0,0,2 
^1,0,1,0 

^1,0,0,1 
^0,1,0,1 





2 
















m 










m 




' 3P 


4 



mr-- 
^1,0,0,1 



+ 



r 
3/ 
mr^ 



/2 2 



m r 



61 



mr^ 



G 



1,0,0,1 



G 



2,0,0,0 



As a simplification we consider an expansion in equations (16)-(21) up to second order 
in quantum variables. The corresponding system is 

(22) 
(23) 
(24) 
(25) 
(26) 
(27) 

(28) 
(29) 



0,2,0,0 



+ 



3/2 
2mr 



:G 



2,0,0,0 



21 



mr^ 



G 



1,0,0,1 



2mr 
1 



1,1,0,0 

- k 

^0,0,0,2 



1,1,0,0 



4/ 



mr^ 



G 



0,1,0,1 



mr^ 
2 

mr^ 



+ 



21 



:G 



0,0,1,1 



+ 



mr 
Al 

mr-^ 



^1,0,0,1 



1,0,1,0 



G 



0,1,1,0 



m 

^0,1,0,1 
m ' 

2mr 
1 



mr^ 



21 



mr-' 



G 



2,0,0,0 



1,0,0,1 



21 



:G 



0,0,0,2 



mr'' 



(30) 
(31) 

(32) 
(33) 



G 



0,1,0,1 



mr^ 



3/2 
2mr 



:G 



1,0,1,0 



+ 



21 



mr-- 



1,1,0,0 _ ^0,0,1,1 
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(34) 

where we already inserted po = 1 = const.. 

One can see that if all the moments G"'''''^''^ are set to zero, we recover the classical 



equations of motion. Equation (30) can immediately be solved G'^''^'^'^ = AZ^ = const. 
Indeed, from (21 ) it can be seen that if ai = 6i = 02 = 0, 62 = then G'^'^'^'^a is always 
zero and all the dispersions of the angular momentum are constants, A/" = const. With 
this simplification the system reduces to a set of twelve coupled differential equations 
for classical and quantum variables. It is evident that we will not find non-spreading 
solutions as in the harmonic oscillator case [21 [10] . 

4. Evolution and quantum back-reaction 

In this section we analyze the evolution and quantum-corrected behavior of the system 



under consideration according to equations (22)- (34). We discuss two cases of interest 



4.1. Zero angular momentum 

The case / = corresponds to the so called one dimensional Hydrogen atom. Although 
one might think that this system as physically uninteresting it actually has interesting 
features [12]. It has been used to model Hydrogen atoms in the presence of strong 
magnetic fields as in astrophysical systems; it is also useful in modeling Rydberg atoms 
in external fields, the behavior of certain electrons near the surface of helium [13] . 
Classically this case reduces to a one-dimensional problem with radial equations: 

Pr . k 

r = — , pr = ^. (35) 

As it evolves in time the position of the particle begins to increase as its velocity 
decreases, after that the radius starts to decrease as the velocity increases in the opposite 
direction (see Fig. [T|. In Fig. [2| we show the phase space of the system in which these 
features are clear. 

However, it turns out that the quantum corrected behavior, up to second order in 
the moments, is not reduced to a one effective dimension when / = 0. As we can see 



from equations (22)-(34), by setting / = 0, and even A/" = 0, there are still quantum 



fluctuations related to the angular momentum as C""'^® and G^' '^*, from equations (32) 



and (33), which are not necessarily zero. This is not surprising since, from (g, we see 
that the term G'^'P* is 

Qr^Pe = Gi'O'O'i = ((f - r){pe - pe))^,^, = (f pe) - Ir, (36) 



and similarly for G^'''^"; these terms do not vanish. Indeed, from equation (24) we 
notice that 6 does not vanish but is proportional to G^''^" . This implies that quantum 
fluctuations induce a two-dimensional motion in the system, modifying the classical 
behavior of r and Pr, as shown in Fig. [T]and[2j The two dimensional quantum-corrected 
behavior of 9{t) is shown in Fig. |3^. In Fig. |3}d it is clear that the actual orbit is 
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(a) (b) 

Figure 1. (a) The position r{t) and (b) the momentum Pr{t) as functions of time, 
for I — 0. The dashed (red) hue corresponds to the classical behavior for m = fc = 1. 
The other curves correspond to the second order effective quantum behavior; when 
A/^ = 0, solid (blue) curve and for H^P ^ dot-dashed (blue) curve. 

Pr 



Figure 2. Classical and quantum corrected phase space for I — 0. The dashed (red) 
is the classical behavior, the quantum evolution for — 0, solid (blue) curve and for 
A/2 ^ dot-dashed (blue) curve. 



not one-dimensional but two-dimensional. This is a purely quantum effect that has no 
classical parallel. 

If we consider / 7^ in the quantum-corrected equations of motion, we notice that 
the system is effectively two dimensional in the classical sense. If we look at equation 
(24) we see that in addition to the term C'^*, there are also two other contributions 
proportional to C'^r"^ and r"^ not present when I = 0. The expectation value 
of products of linear and angular momentum operators in general do not commute, 
therefore the angular equation does not vanish. We see here that, at the semiclassical 
effective level, the classical observables acquire quantum corrections that are not 
negligible and that back-react latter on the classical variables. 

This system does not correspond exactly to the effective description of the one 
dimensional Hydrogen atom. To be so all moments related to angular variables would 
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necessarily equate to zero. Amusingly, when we impose this in our analysis we notice that 
the corresponding solution has qualitatively the same behavior as in Fig. [1} although 
in this case the motion is purely one dimensional. 

4-2. Full two dimensional case 

In standard quantum mechanics energy levels are obtained by solving the radial 
Schrodinger equation and imposing asymptotic conditions on the solution [TT]. 
Observables Q are then obtained from the wave function il){v) (for stationary solutions) 
in the usual way as its expectation values Q = (Q) = J 'ip* Q'ip d^x. 

The approach based in effective equations treats the dispersions as quantum 
corrections to the evolution. Because dispersions also evolve in time, we must analyze 
whether our considerations are valid, i.e., whether the quantum variables can be 
considered as perturbations, and under what circumstances. Our main goal is to study 
the evolution of expectation values and its deviation from the classical picture, then we 



now proceed to solve numerically equations (22)- (34). 

In Fig. |§i it can be seen that for initial conditions that consider small enough 
dispersions, there is a region close to t = where the dispersions are small compared 
to the classical variables, allowing a perturbative evolution. Furthermore, we notice 
in Fig. |4]d, that for large times the radial coordinate oscillates with small amplitude 
around r = 1. The angular variable and its dispersion have a similar behavior. 

Note that, as the fluctuations evolve, they begin to grow until they are comparable 
or even bigger to the expectation values. In this region they could no longer be 
considered as perturbations. 

For Pr/^ G'P'-P'- we can see a similar oscillatory evolution in Fig. sj for the same 
initial conditions. In the first part a perturbative approximation can be implemented 
while for larger times, it oscillates around zero with small amplitude. The effective radial 
momentum as compared with the classical one, starts close to it but as time increases 
the momentum becomes larger as well as its oscillation frequency. Again this implies 
that for large times the dispersion are of the same magnitude as the momentum. 
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(a) 



(b) 



Figure 4. (a) Plot of r jyfG^ with initial conditions that consider small enough values 
for dispersions, (b) The radial coordinate as a function of time. Compared with the 
classical behavior, the effective one has a decreasing amplitude as evolves, and for large 
times oscillates around r = \. In this plots we consider m = 1, fc = 2 and I = \. The 
initial conditions for unperturbed curves are r = 1, pr — 1 and 9 — t:; and for the 
moments we consider all of order 0.01. 




(a) 



(b) 



Figure 5. (a) Plot of Pr/\/ Gpp and (b) of the radial momentum, both with the same 
initial conditions as in Fig. |4j 



One can always impose initial conditions complying with the uncertainty principle 
(|2| for each pair of conjugate variables. We can see this in Fig. [6] for the radius, Fig. [7] 
for the ratio of the corresponding dispersions and Fig. [8] for their momenta, taking 
for the dispersions initial conditions progressively smaller. Although the behavior is 
qualitatively similar to that shown in Fig. |4] and Fig. |5| we can notice that, for the 
initial conditions that fulfill the uncertainty relations, we have better control of the 
approximation. This is because the effective evolution is very close to the classical one. 
Indeed, the smaller the initial condition for the perturbation, the longer the effective 
trajectory remains close to the classical one. One can see this kind of behavior, where the 
classical orbits decay after some time, when one considers the hydrogen atom interacting 
with external fields [17j. 

It is expected that the energy of the system is also back-reacted, which indeed is 



verified from the definition of the effective quantum energy (16). However, as classically 
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Figure 6. Plots of the radius as a function of time for initial conditions fulfilling 
the uncertainty relation (a) Dispersions of order 0.001, (b) order 10""* and (c) 
10~^. We can see that the smaller the perturbations, the longer they remain near the 
classical trajectory, allowing longer perturbation regions. 




Figure 7. Plots of the ratio r/vG^ as a function of time for the same initial conditions 
as in Fig. [6] Interestingly, for these initial conditions the graph presents some peaks 
which indicates that dispersions may be considered as perturbations around those 
regions. 



P 




(a) (b) 



Figure 8. (a) Plot of pr and (b) Pr/^ Gpp, for initial conditions of dispersions of order 
10^^. There are piecewise regions where the perturbative approximation is well suited. 



the energy is a constant of motion, in the perturbative regime the effective behavior 
of the corrected energy oscillates about the classical energy with small amplitudes that 
depend on the initial conditions for the dispersions, as shown in Fig. [9j 
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Figure 9. The classical (dashed red line) and quantum corrected (solid blue line) 
energy as function of time. The initial conditions for dispersions are of order 0.001. 



Let us recall that classically the solution of the Kepler problem corresponds to a 
conic-section orbit whose eccentricity and perihelion depends on the energy and angular 
momentum. Although in the quantum-corrected study that we are discussing (as in any 
quantum/statistical description) one cannot have a closed analytical expression for the 



effective orbit, one can solve the system (22)-(34) numerically. With initial conditions 



fulfilling the uncertainty relations, we obtain an almost elliptical (open) orbit that is 



modified by the effect of back-reaction of quantum variables, as shown in Fig. 10 Under 
these conditions the effective orbit starts very close to the classical one and, as it evolves, 
quantum effects increasingly drive it away. 





(a) 



(b) 



Figure 10. The classical (dashed red line) and quantum corrected (solid blue line) 
orbits for dispersions of order (a) 10""* and (c) 10^^. Classical and quantum-corrected 
orbits starts close, but at large times the effective one spreads away. 



The departure of the effective orbit from the classical one can also be understood 
by looking at the uncertainties in r. These are the curves located at each side of 
the quantum trajectory r(6'), (i.e. the brown and purple thin curves), in Fig. 11 



Here we notice how, at the beginning of evolution, both the effective and classical 
orbits are within the range of dispersions ±-\/G'''' = ±Ar. As the system evolves 
we can see how the uncertainties start to grow and disperse. Interestingly at some 
point the uncertainties decrease and we recover the effective trajectory. For large initial 
conditions of the uncertainties this does not happen very close to the classical trajectory, 
but as one decreases the initial values this happens closer to the classical behavior. 
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Fig. 11 Actually, from the previous discussion we know that there is a region where the 
magnitude of the dispersion is comparable to that of the expectation values. Then we 
can interpret these as the quantum effects that keep the orbit away from the classical 
trajectory for large times. This behavior is similar to the one in [18j, where minimum 
uncertainty states were constructed in which the wave function can be described by 
classical equations for short times, while for longer times the wave packet that lives in 
the elliptical trajectory begins to spread. 




(a) (b) (c) 

Figure 11. Plot of the quantum-corrected (blue) orbit surrounded by the 
corresponding dynamical uncertainties r ± Ar (thin curves), for initial conditions of 
the dispersions of order (a) 0.001, (b) 10~^ and (c) 10~^. Also the classical orbit (red 
dashed curve) is shown. Initially the uncertainties are close to, and at each side of, 
the effective orbit, then they start to disperse. At large times they diverge completely 
from the classical behavior. 



We can further appreciate the quantum evolution of the system by looking at 
the corresponding radial phase space diagram. Classically, for a particle rotating in 
an elliptic orbit in configuration space, the radial phase diagram is a closed curve 
representing a bound state in which the particle attains a maximum and a minimum 
value of its radius from the force center (dashed (red) curve in Fig. 125,). The solid 
curve in Fig. 12 a, corresponds to the quantum-corrected evolution; we notice that the 
effective quantum behavior starts close to the classical diagram and then disperses. 
This diagram shows how the quantum effects make the effective system depart from the 
classical one. As time evolves we see that the radius becomes increasingly smaller and 
the radial momentum grows. Moreover, after a long enough period of time, we can see 
in Fig. 12 D that, as the radius gets locahzed, the radial momentum disperses: this is a 
direct manifestation of the uncertainty principle in phase space. 

Fig. 13 shows how the dispersions evolve and to what extent the uncertainty 
principle is fulfilled. From Fig. 13 we see that the inequality is satisfied for short times, 
during the phase where the expansion can be treated perturbatively, as mentioned above. 
There is only one region where the product presents a peak but then remains flat for 
some more time. For large times the perturbative assumption is no longer valid, as 
the product of uncertainties grow larger and larger. This is a similar behavior as the 
one presented by a radial squeezed state studied in [Tl|, where the uncertainty product 
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(a) (b) 

Figure 12. (a) The dashed (red) hne corresponds to the classical radial phase diagram 
while the solid (blue) curve is the effective quantum diagram for initial conditions of 
dispersions of order 0.01. (b) For later times the radial coordinate is fully localized 
while the radial momentum disperse completely. 




features a cyclic behavior. In the same way we can say that our system begins in a state 
with minimum uncertainty, not squeezed, and for larger times it gets dispersed. 

Ap Ar 




(a) (b) 

Figure 13. The uncertainty relations (a) ArAp,. and (b) G"''' G'P'-P'- - {C^'-f as 
function of time (order 10^^). During the perturbative stage, the uncertainty relations 
are satisfied; for larger times the inequality is violated and quantum effects become 
increasingly dominant. 



4-. 3. Different approaches to semi- classical limit of the Hydrogen atom 

Even though the non-relativistic Hydrogen atom has been solved exactly it was latter 
proved that the WKB approximation was rather poor in this case, which is usually 
attributed to the singular behavior of the potential at the origin r = 0. It was shown in 
[I5] that WKB can be performed in r — > with precise results if the centrifugal term 

2mr^ 

is replaced by 



Vr 



c 



2mr^ 
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which in turn gives very accurate energy eigenvalues for the WKB approximation, even 
at lowest order. The main difference between WKB and the method presented here 
is the scope of the treatment: while both methods intend to establish a semiclassical 
approach to the quantum mechanical system at hand, the WKB aims at obtaining, by 
means of an series expansion on the wave function and corrections to the energy 
levels. That is, one writes the wave function as [1] 



and the WKB approximation consists in expanding W{y) in powers of h and in neglecting 
in the Schrodinger equation terms of order or higher. Besides the fact that the method 
aims in approximating the wave function it fails in the vicinity of the classical turning 
points, i.e., E = V{x). To circumvent this one defines the connection functions in order 
to connect the wave functions in the regions adjacent to the turning points. 

Evolution and description of the quantum mechanical system in our method is not 
based on a wave function but rather in its momentous description and their quantum 
back-reaction on expectation values. One can compute the energy of the system and, 
probably the most interesting aspect of our description, also the semiclassical evolution 
of observables and orbits at different regimes in the parameter space. However, one 
feature that both schemes (and any approximation method in quantum mechanics) 
share is the growing in time of quantum corrections that render the approach invalid as 
well as the oscillatory character of the evolution in our case or in the expression of wave 
function and energy in the WKB. This is a well known effect in quantum mechanics. 

The effective equations method presented here is an alternative method that in 
principle should reproduce all the results of standard quantum mechanics as long as 
one consider all the infinite moments. As mentioned, this method is quite suitable for 
a semi-classical analysis, defining a region where the moments can be considered as 
perturbations. 

The semi-classical states of the Hydrogen atom have been explored in terms of the 
correspondence principle for large principal quantum number n. In such a case a wave 
packet can follow a classical orbit. Interestingly, given the conservation of energy and 
angular momentum, there are no transversal spreading in such states. The simplest 
of this states has a dispersion of order and it evolves in such a way that the 

wave packet is distributed almost uniformly on a circular orbit [16j. The dynamics of 
these states have both classical and quantum features: for short times the motion is 
classical, after that the quantum dynamics becomes dominant [T^. Something very 
similar happens in our case, although we do not necessarily describe the same state, as 
the effective equations method provides us with minimal uncertainty dynamical states. 
There is also a formalism to construct radial and angular coherent and squeezed states 
that show the main features of classical motion p3l [19] . Semi-classical states have also 
been constructed from the coherent states of a 4 dimensional oscillator which is reached 
by the KS mapping of the 3 dimensional Hydrogen atom pO] . 
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It is interesting to study the semi-classical states of the hydrogen atom as it has 
been possible to find them experimentally in the laboratory through the excitation of 
the outermost electrons of certain atoms, with ultrashort laser pulses [21]. Recently, 
there have been obtained experimentally new localized wave packets in very large n 
states that travel in nearly circular orbits [22] . 

The analysis presented here provides an alternative tool for the study of some 
different semi-classical states of the Hydrogen atom that had not been considered before, 
indeed we can see that our states are not of minimum uncertainty for all times. This 



is even clearer if we consider the ratio Ar/Ap^ in Fig. 14 One can see that there 
is a group of peaks indicating that at those times the uncertainty in r increases as 
uncertainty in pr decreases, while in the flat regions the opposite occurs. As in Fig. 13 
we can compare this with the behavior of the radial squeezed state [H] , where this ratio 
starts with a large amplitude and then damped oscillates until it reaches a minimum 
although continues oscillating, showing the squeezing in p^.- Clearly our state is not a 
squeezed one, but is the closest one that begins satisfying uncertainty relations. 

Ar 
Ap 



t 



Figure 14. Ratio Ar/Ap^- The peaks indicate that there is more spread in the radius 
than in the momentum at that time, i.e. states are squeezed for short periods. 



5. Discussion 

In this work we have obtained an effective description of the Hydrogen atom by means of 
a momentous formulation of quantum mechanics. We were able to obtain a semiclassical 
picture of the evolution of the classical dynamical variables by studying their equations 
of motion, back reacted by quantum dispersions and spreadings. This is done in regions 
where the moments can be considered as perturbations, which also allows us to make a 
consistent truncation of the equations of motion. 

The equations of motion for expectation values are obtained from an effective 
Hamiltonian that acquires quantum corrections represented by spreadings to all orders. 
We have seen that the system reduces to the Kepler problem when one sets all the 
quantum fluctuations to zero. The case of zero angular orbital momentum, / = 0, 
that classically corresponds to a one dimensional system, gets so strongly quantum 
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corrected that it effectively corresponds to a two dimensional system. This behavior 
has no classical analogous, and in this way is a similar effect as the quantum tunneling 
phenomena. 

The evolution of the system is controlled by generalizations of the Heisenberg 
uncertainty relations among quantum variables that are satisfied during the perturbative 
period, thus providing us with a consistent analysis and evolution. 

As we pointed out we were able to determine the evolution of expectation values 
of variables which enabled us to determine an effective orbit for the electron, which is 
close to the classical elliptical orbit for short initial times. For larger times it is diverted 



to an open quasi-elliptical orbit as shown in Figs. 10, 11 and 12 One can see from 
these figures how quantum effects modify the classical behavior, realizing that, as the 
system evolves, there is a switching between the radial and momentum spreadings. It 
is interesting to point out that the behavior of the orbits and the spreadings shown in 
Fig. [TT| although not identical, are similar to those reported in the literature where the 
wave packets follow elliptical orbits for a while and then spread. 

Although this method gives a semi-classical behavior in a well-defined region, the 
states that we use are not the minimum uncertainty states that are usually called 
coherent states for the Hydrogen atom, or the Rydberg-type states similar to those 
that have been experimentally produced through ultra-short laser pulses. Our states 
are new states describing the semiclassical behavior of the Hydrogen atom. Thus we 
emphasize the qualitative features of this analysis, the method of effective equations 
and its robustness in obtaining new states suited to describe the semiclassical behavior. 

Even though we were able to solve numerically the equations of motion, and with 
this obtaining the evolution of the system providing a rather interesting description, it 
is still needed to obtain a set of physical initial conditions that will allow us to compare 
quantitatively our results with experimental or theoretical data in a consistent manner. 
Due to the very complicated nature of the dynamical system and the infinite number of 
quantum variables this task is yet under further study. 

Due to the complexity of the effective variables and their equations of motion we 
restricted ourselves to the two dimensional case, corresponding to the classical version, 
and then it is necessary to explore the three dimensional system. Given the results 
obtained in the one dimensional case, where quantum effects induce a two dimensional 
behavior, we could obtain modifications to the two dimensional motion due to quantum 
back-reaction by considering the remaining angular variable and its associated momenta, 
fixed in the present analysis. In such case the corresponding results should be more 
accurate and closer to those currently reported in the literature and experimentally. 

The implications and possible applications of our analysis to other systems and 
phenomena are promising. For instance, the implementation of this method to 
Hydrogen-like atoms and systems should be straightforward and may result in a more 
detailed description and a deeper understanding of those systems. Note that, even for 
the Helium atom (a system with two electrons), there is no exact solution, and any 
study of this system is performed at an approximation level, being ours one that could 
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Appendix A. Higher order moments 

As shown in fS] , as one considers higher order moments the average number of terms in 
the equations increases exponentially. Here we show the system of effective evolution 
equations up to third order, that corresponds to thirty coupled equations 
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We omit equations for r, pe, G^'^''^'^ and G^'°'°'^ because they are not modified with 
respect to previous expressions of second order. Also we replace G°'°'°'2 with AP. 
Beside this system there are another eighteen equations: 
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From (|A.23|) we immediately see that G"'*^'"'^ is a constant, as we already noticed 



for the moments of G^'^'^'"" to any order n. 
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